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Numerical simulations of the incompressible Euler equations are performed using the Taylor- 
Green vortex initial conditions and resolutions up to 4096^. The results are analyzed in terms of 
the classical analyticity strip method and Beale, Kato and Majda (BKM) theorem. A well-resolved 
acceleration of the time-decay of the width of the analyticity strip S{t) is observed at the highest 
resolution for 3.7 < t < 3.85 while preliminary 3-D visualizations show the collision of vortex sheets. 
The BKM criterium on the power-law growth of supremum of the vorticity, applied on the same 
time-interval, is not inconsistent with the occurrence of a singularity around t ~ 4. 

These new findings lead us to investigate how fast the analyticity strip width needs to decrease to 
zero in order to sustain a finite-time singularity consistent with the BKM theorem. A new simple 
bound of the supremum norm of vorticity in terms of the energy spectrum is introduced and used to 
combine the BKM theorem with the analyticity-strip method. It is shown that a finite-time blowup 
can exist only if S{t) vanishes sufficiently fast at the singularity time. In particular, if a power law is 
assumed for S{t) then its exponent must be greater than some critical value, thus providing a new 
test that is applied to our 4096^ Taylor-Green numerical simulation. 

Our main conclusion is that the numerical results are not inconsistent with a singularity but that 
higher-resolution studies are needed to extend the time-interval on which a well-resolved power-law 
behavior of 5{t) takes place, and check whether the new regime is genuine and not simply a crossover 
to a faster exponential decay. 

PACS numbers: 47.10.A,47.11.Kb,47.15.ki 



I. INTRODUCTION 

A central open question in classical fluid dynamics 
is whether the incompressible three-dimensional Euler 
equations with smooth initial conditions develop a singu- 
larity after a finite time. A key result was established in 
the late eighties by Beale, Kato and Majda (BKM). The 
BKM theorem |T] states that blowup (if it takes place) 
requires the time-integral of the supremum of the vortic- 
ity to become infinite (see the review by Bardos and Titi 
[2]). Many studies have been performed using the BKM 
result to monitor the growth of the vorticity supremum 
in numerical simulations in order to conclude yes or no 
regarding the question of whether a finite-time singular- 
ity might develop. The answer is somewhat mixed, see 
e.g. references [5HS] Eind the recent review by Gibbon 
[6]. Other conditional theoretical results, going beyond 
the BKM theorem, were obtained in a pioneering paper 
by Constantin, Fefferman and Majda 0. They showed 
that the evolution of the direction of vorticity posed ge- 
ometric constraints on potentially singular solutions for 
the 3D Euler equation [7]. This point of view was further 
developed by Deng, Hou and Yu in references [5] and [p]. 

An alternative way to extract insights on the singular- 
ity problem from numerical simulations is the so-called 
analyticity strip method [10]. In this method the time 
is considered as a real variable and the space-coordinates 



are considered as complex variables. The so-called "width 
of the analyticity strip" 6{> 0) is defined as the imaginary 
part of the complex-space singularity of the velocity field 
nearest to the real space. The idea is to monitor S{t) as a 
function of time t. This method uses the rigorous result 
[11] that a real-space singularity of the Euler equations 
occurring at time T* must be preceded by a non-zero 6{t) 
that vanishes at T*. Using spectral methods [T^, 5{t) is 
obtained directly from the high-wavenumber exponential 
fall off of the spatial Fourier transform of the solution 
. This method effectively provides a "distance to the 
singularity" given by 6{t) |14| . which cannot be obtained 
from the general BKM theorem. 

Note that the BKM theorem is more robust than the 
analyticity-strip method in the sense that it applies to 
velocity fields that do not need to be analytic. How- 
ever, in the present paper we will concentrate on initial 
conditions that are analytic. In this case, there is a well- 
known result that states: In three dimensions with peri- 
odic boundary conditions and analytic initial conditions, 
analyticity is preserved as long as the velocity is contin- 
uously differentiable (C^) in the real domain [11_. The 
BKM theorem allows for a strengthening of this result: 
analyticity is actually preserved as long as the vorticity 
is finite [14 . 

The analyticity-strip method has been applied to probe 
the Euler singularity problem using a standard periodic 
(and analytical) initial data: the so-called Taylor-Green 
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(TG) Vortex [TS' . We now give a short review of what is 
already known about the TG dynamics. Numerical sim- 
ulations of the TG flow were performed with resolution 
increasing over the years, as more computing power be- 
came available. It was found that except for very short 
times and for as long as S{t) can be reliably measured, 
it displays almost perfect exponential decrease. Simula- 
tions performed in 1982 on a grid of 256^ points obtained 
S{t) - 2.60e-*/"-57 (for t up to 2.5) [H]. This behavior 
was confirmed in 1992 at resolution 864^ [17]. More than 
20 years after the first study, simulations performed on a 
grid of 2048^ points yielded S{t) - 2.70e~*/°-^6 (for t up 
to 3.7) ^18J. If these results could be safely extrapolated 
to later times then the Taylor-Green vortex would never 
develop a real singularity [13' . 

The present paper has two main goals. One is to re- 
port on and analyze new simulations of the TG vortex 
that are performed at resolution 4096^ . These new simu- 
lations show, for the first time, a well-resolved change of 
regime, leading to a faster decay of d{t) happening at a 
time where preliminary 3-D visualizations show the colli- 
sion of vortex sheets [26 . The second goal of this paper 
is to answer the following question, motivated by the new 
behavior of the TG vortex: how fast does the analyticity- 
strip width have to decrease to zero in order to sustain 
a finite-time singularity, consistent with the BKM theo- 
rem? To the best of our knowledge, this question has not 
been formulated previously. 

To answer this question we introduce a new bound 
of the supremum norm of vorticity in terms of the en- 
ergy spectrum. We then use this bound to combine the 
BKM theorem with the analyticity-strip method. This 
new bound is sharper than usual bounds. We show that 
a finite-time blowup exists only if the analyticity-strip 
width goes to zero sufficiently fast at the singularity time. 
If a power-law behavior is assumed for S(t) then its expo- 
nent must be greater than some critical value. In other 
words, we provide a powerful test that can potentially 
rule out the existence of a finite-time singularity in a 
given numerical solution of Euler equations. We apply 
this test to the data from the latest 4096'^ Taylor-Green 
numerical simulation in order to see if the change of be- 
havior in S{t) can be consistent with a singularity. 

The paper is organized as follows: Section |TT] is de- 
voted to the basic definitions, symmetries and numerical 
method related to the inviscid Taylor-Green vortex. 

In Sec. |III| the new high-resolution Taylor-Green re- 
sults are presented and axe analyzed classically in terms 
of analyticity-strip method and BKM. 

In Sec. |IV[ the analyticity-strip method and BKM 
Theorem are bridged together. The section starts with 
heuristic arguments that are next formalized in a math- 
ematical framework of definitions, hypotheses and theo- 
rems. 

In Sec. |V] our new theoretical results are used to an- 
alyze the behavior of the decrement. 
Section IVTl is our conclusion. 

The generalization to non TG-symmetric periodic 



fiows of the results presented in Sec. |IV| are described 
in an appendix. 

II. DEFINITION OF THE SYSTEM 

A. Basic definitions 

Let us consider the 3D incompressible Euler equa- 
tions for the velocity field u{x,y,z,t) G defined for 
{x, y, z) g and in a time interval t e [0, T): 

du , , 

— +u-Vu = -Vp, V-u = 0. (1) 

The Taylor-Green (TG) flow [H] is defined by the 2tt- 
periodic initial data u(a;, y, z, 0) = u'^^(a;, y, z), where 

u'^^ — (sin(a;) cos{y) cos(z), — cos(a;) sm{y) cos(z), 0). 

The periodicity of u allows us to define the (standard) 
Fourier representation 

u(k, t) = / u(x,t) exp(-ikx)d^a; (2) 

(27r) Jjj 

u(x,t) = ^ u(k, i) exp(ikx), (3) 

The kinetic energy spectrum E{k,t) is defined as the 
sum over spherical shells 

mt)^l J2 is(k,oP, (4) 

fc-l/2<|k|<*: + l/2 

and the total energy 

''^2(i)^/j"^^'^^''''^^2^'"^''^^''' 

is independent of time because u satisfies the 3D Euler 
equations {[TJ. 

B. Symmetries 

A number of the symmetries of u"""^ are compatible 
with the equation of motions. They are, first, rotational 
symmetries of angle tt around the axis (x = z = 7r/2) 
and [x = z = 7r/2); and of angle 7r/2 around the axis 
{x = y — 7t/2). a second set of symmetries corresponds 
to planes of mirror symmetry: a; = 0, tt, y = 0, tt and 
z = 0, TT. On the symmetry planes, the velocity u"^*^ and 
the vorticity lj^'^ = V x u'^^ are (respectively) parallel 
and perpendicular to these planes that form the sides of 
the so-called impermeable box which confines the fiow. 

It is demonstrated in reference [16 that these sym- 
metries imply that the Fourier expansions coefficients of 
the velocity field in Eq. ([3| u{m,n,p,t) vanishes unless 
m,n,p are either all even or all odd integers. This fact 
can be used in a standard way [16] to reduce memory 
storage and speed up computations. 
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C. Numerical method 

The Euler Equations ([T]) are solved numerically us- 
ing standard [T2] pseudo-spectral methods with resolu- 
tion N. Time marching is done with a second-order 
Runge-Kutta finite-difference scheme. The solutions are 
dealiased by suppressing, at each time step, the modes for 
which at least one wave-vector component exceeds two- 
thirds of the maximum wave- number N/2 (thus a 4096^ 
run is truncated at fc > /cmax = 1365). 

The simulations reported in this paper were performed 
using a special purpose symmetric parallel code devel- 
oped from that described in [TQ] [20] . The workload for 
a timestep is (roughly) twice that of a general periodic 
code running at a quarter of the resolution. Specifically, 
at a given computational cost, the ratio of the largest 
to the smallest scale available to a computation with en- 
forced Taylor-Green symmetries is enhanced by a factor 
of 4 in linear resolution. This leads to a factor of 32 
savings in total computational time and memory usage. 
The code is based on FFTW and a hybrid MPI-OpenMP 
scheme derived from that described in |j21j. The runs 
were performed on the IDRIS BlueGene/P machine. At 
resolution 4096'^ we used 512 MPI processes, each process 
spawning 4 OpenMP threads. 



III. NUMERICAL RESULTS AND CLASSICAL 
ANALYSIS 

A. Energy spectra, maximum vorticity and 
collision of vortex sheets 
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FIG. 1: (Color online) Temporal evolution of TG flow: 
a) energy spectra E(k,t) (see Eq. Q) at t = 
(1.3,1.9,2.5,2.9,3.4,4.0) and b) maximum of vorticity 
^)lloo- Results from runs performed at different reso- 
lutions are displayed together: 512^ (brown triangles), 1024^ 
(blue squares), 2048^ ( green diamonds) and 4096"^ (red cir- 
cles). 



Runs were performed at resolutions 512^^, 1024^, 2048'^ 
and 4096^ 

The behavior of the energy spectra Q and the spatial 
maximum of the norm of the vorticity a; — V x u are 
presented in Fig. [l] 

It is apparent in Fig. [ija) that resolution-dependent 
even-odd oscillations are present, at certain times, on 
the TG energy spectrum. Note that this behavior is 
produced when the tail of the spectrum rises above the 
round-off error ~ 10"'^^. This phenomenon can be ex- 
plained in terms of a resonance |22| . along the lines devel- 
oped in reference [23| . In practice we will deal with this 
problem by averaging the spectrum over shells of width 
Afc = 2. Apart from this it can be seen that spectra com- 
puted using different resolutions are in good agreement 
for all times. 

In contrast, it is visible in Fig. [ijb) that the maxi- 
mum of vorticity ||'^(', OIloo computed at different res- 
olutions are in agreement only up to some resolution- 
dependent time (see the inset). The fact that ||a;(-,t)||^ 
at a given time t > 3.7 decreases if one truncates the 
higher wavenumbers of the velocity field (see Fig. [ijb)) 
strongly suggests that ||a;(-, i)|j^ has significant contribu- 
tions coming from high- wavenumbers modes. This forms 



the basis of the heuristic argument presented below in 

Sec. nvn 

Figure [2] shows 3D visualizations (using the VAPOR 
|27j software) of the high vorticity regions in the imper- 
meable box, corresponding to the 4096'^ run at late times. 
A thin vortex sheet is apparent in Figj2ja) on the vertical 
faces X = 0, TT and y = 0, tt of the impermeable box. 

The emergence of this thin vortex sheet is well under- 
stood by simple dynamical arguments about the fiow on 
the faces of the impermeable box that were first given in 
reference [16]. We now briefly review these arguments. 
The initial vortex on the bottom face is first forced by 
centrifugal action to spiral outwards toward the edges 
and then up the side faces. A corresponding outfiow on 
the top face and downflow from the top edges onto the 
side faces leads to a convergence of fluid near the hori- 
zontal centreline of each side face, from where it is forced 
back into the centre of the box and subsequently back 
to the top and bottom faces. The vorticity on the side 
faces is efficiently produced in the zone of convergence, 
and builds up rapidly into a vortex sheet (see Figs. 1 and 
2 of reference [H] and Fig. 8 of reference [17]). 

While these considerations explain the presence of the 
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FIG. 2: (Color online) 3D visualization of TG vorticity | V x u| 
at resolution 4096'': a) full impermeable box < 2: < tt, 
< y < 71" and 0<z<-Katt — 3.75. Zooms over the 
subbox marked near x — y = n, z — tt/2 are displayed in b) 
a.t t = 3.5, in c) a.t t — 3.75 and in d) at t = 4.0. 



thin vortex sheet in Figj2ja), the dynamics presented in 
Figj2jb-d) also involves the collision of vortex sheets hap- 
pening near the edge x = y = tt, close to z = tt/2. Note 
that, as stated above in Sec. |II B| the vortex lines are per- 
pendicular to the faces of the impermeable box. Thus, 
because the collision takes place near an edge, the corre- 
sponding vortex lines must be highly curved, with strong 
variations of the direction of vorticity. The geometric 
constraints on potential singularities posed by the evolu- 
tion of the direction of vorticity developed in references 
[THS] could be applied to the situation described in Fig. 
[|] Howe ver, such an analysis goes beyond the BKM the- 
orem and involves extensive post-processing of very large 
datasets. This task is thus left for further work and we 
concentrate here on simple BKM diagnostics for the vor- 
ticity supremum and analyticity strip analysis of energy 
spectra. 



B. Analyticity-strip analysis of energy spectra 

The analyticity-strip method [TD] is based on the fact 
that when the velocity field is analytic in space the energy 
spectrum satisfies E{k, t) the asymptotic 

'ultraviolet region' A: 3> 1, with a proportionality factor 
that may contain an algebraic decay in fc, a multiplicative 



function of time and, depending on the complexity of the 
physical fiow, even an oscillatory (in k) modulation |lS] . 
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FIG. 3: (Color online) Comparison of fit ([5| (solid black line) 
and spectrum at resolution 4096^ (red markers); times and fit 
intervals are indicated in the legend. 



The basic idea is thus to assume that E{k,t) can be 
well approximated by a function of the form 

E{k,t) « C(i)fc-"(*) e-^kSit) ^ 

in some wave numbers interval between 1 and fcmax = 
[A'^/3J (the maximum wavenumber permitted by the nu- 
merical resolution A'^). The common procedure to de- 
termine C{t),n{t),S{t) is to perform a least-square fit 
at each time t on the logarithm of the energy spectrum 
E{k,t), using the functional form 

\nE{k, t) = In C(<) - n{t) In k-2kS{t) . (5) 

The error on the fit interval ki < k < k2, 

X^{t) = {\nE{k,t) -\nClt) + n{t) \nk + 2k6{t)f , 
k—ki 

is minimized by solving the equations dx^/dC = 0, 
dx^/dn = and dx^/dS = 0. Note that these equations 
are linear in the parameters InC(i), n{t) and S{t). 

The transient oscillations of the energy spectrum ob- 
served at the highest wavenumbers (see above Fig. [l|^a) 
are eliminated by averaging the TG spectrum on shells 
of width Afc — 2 before performing the fit ^16j . 

We present in Fig. [3] examples of TG energy spec- 
tra fitted in such a way on the intervals 2 < k < 
min(A:*, fcmax), where k* — uiiE(k)<io-^^{k) denotes the 
beginning of round off noise. It is apparent that the fits 
are globally of a good quality. 
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FIG. 4: (Color online) Time evolution of energy spectrum fit parameters (see Eq. (jSjl and Fig. |3[a)) constant C; b) prefactor 
n; c) decrement 5 (horizontal lines : ^fcmax ~ 2, dashed black line: exponential law a) decay rate —d{ln{5{t))/dt. 

Results corresponding to different resolutions are displayed together: 5f2^ (brown triangles), f024'' (blue squares), 2048^ 
(green diamonds) and 4096^ (red circles). 



The time evolution of the fit parameters C, S and n 
computed at different resolutions are displayed in Fig. [4] 
The measure of the fit parameters is reliable as long as 
d(t) remains larger than a few mesh sizes, a condition 
required for the smallest scales to be accurately resolved 
and spectral convergence ensured. Thus the dimension- 
less quantity ^fc^ax is a measure of spectral convergence. 

It is conventional [TB] to define a 'reliability time' Trci 
by the condition 

'5(Trcl)fcmax = 2, (6) 

and to say that the numerical simulation is reliable for 
times t < Tiei. This reliability time can be extended 
only by increasing the spatial resolution available for the 
simulation, so the more computer power is available the 
larger is the reliability time. 

The resolution-dependent reliability condition ([6| is 
marked by the horizontal lines in Fig. |4|c). The ex- 
ponential law 

,5(t) - 2.70e-*/°-5^ (7) 



that was previously reported at resolution 2048"^ in refer- 
ence [13 is also indicated in Fig. Qc) by a dashed black 
line. It is thus apparent that our lower-resolution re- 
sults well reproduce the previous computations that were 
discussed above in Sec. |l] (see text preceding references 

UMi). 

In Table |l] the reliability time ^ obtained from the 
fit parameter 6 of Fig. |4] is compared with the relia- 
bility time stemming from the exponential behavior (|7|. 
It is apparent by inspection of the Table that the relia- 



Resolution 


Trci (exponential law) 


r,ei (fit) 


5f2^ 


3.05 


3.05 


1024^ 


3.43 


3.44 


2048^ 


3.82 


3.75 


4096^ 


4.21 


3.85 



TABLE I: Reliability time Q deduced from the exponential 
behavior Q compared with the reliability time obtained from 
the fit parameter 5 of Fig. [4] 

bility time of our new 4096'^ results is markedly smaller 
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than that deduced from the exponential law ([t]) ; the lat- 
ter wrongly predicts that simulations at this resolution 
should be reliable until t = 4.21. The departure from the 
exponential behavior is also visible on the inset in Fig. 
ic). 

In order to capture this change of behavior more quan- 
titatively the logarithmic decay rate —d\n{S)/dt, com- 
puted using finite differences in time, is displayed in Fig. 
^d). A clear change in trend is apparent around t — 3.7. 
where the logarithmic decay rate abruptly changes from 
a value near 2 to a value near 8. Note that this change of 
behavior happens at a time that is reliable at resolution 
4096'^ (see insets in Fig. ^c) andjijd). Interestingly, this 
time is close to the reliability time of the 2048^ simula- 
tion. Therefore, the new behavior of accelerated decay 
for times t > 3.7 can only be suggested by the 2048"^ data 
and is here demonstrated for the first time by our 4096"^- 
resolution data. This acceleration of the decay rate of 
S{t) is important because if ([T]) could be safely extrapo- 
lated to later times then the Taylor-Green vortex would 
never develop a real singularity [13 . 

Let us conclude this section by showing that the new 
behavior does not depend on the wavenumber interval 
chosen to perform the fits. 

Indeed, by close inspection of the top curve in Fig. [3] 
one can see that a small amount of systematic errors are 
present at the lowest {k < 100) wavenumbers for large 
times. Excluding the lowest wavenumbers from the fits 
results in less errors (data not shown) . In Table |ll] the 
result of fits performed on the subinterval 103 < k < /c,nax 
are compared with those on the full interval 3 < /c < 
^max that was used up to now . It can be checked on 



Time 


n 


n 


10^ X 5 


10^ X 5 




3 ^max 


103 ^max 


3 ^max 


103 ^max 


3.6 


4.07 


3.95 


4.13 


4.22 


3.65 


4.09 


4.05 


3.73 


3.75 


3.7 


4.09 


4.14 


3.36 


3.31 


3.75 


4.09 


4.19 


2.85 


2.76 


3.8 


4.12 


4.29 


2.10 


1.95 


3.85 


4.13 


4.34 


1.41 


1.22 


3.9 


4.09 


4.34 


0.94 


0.71 



TABLE II: Time evolution of fit parameters n and 5 (see 
Eq. ([5|) on full interval 3 < fc < fcmax (same as in Fig. [4| 
compared with fits on subinterval 103 < k < femax- 

the table that the departure from the exponential law is 
not dependent on the interval chosen to perform the fit. 
The values of n are also in agreement with previously 
pubHshed data [18] . 



C. BKM analysis of vorticity maximum 

In this section we look for eventual singular behavior 
by focusing on the time-dependence of the TG data for 
the vorticity suprcmum ||Li;||oo(i) that is displayed above 
in Fig. ^h). The BKM theorem pi states that blowup 



(if it takes place) requires the time-integral of the supre- 
mum of the vorticity to become infinite. Our analysis 
method, first introduced in amounts to look at evi- 
dence of power-law behavior in the numerical time series 
for 1 1 a; I loo (i) and see if the computed exponent is com- 
patible with blowup of the time integral of ||u;||oo(i)- We 
now proceed to briefly recall the method. 

Let f{t) be the quantity to be studied. In order to 
test if it might blow up or go to zero in a finite time, we 
produce, locally in time, fits of power law behavior of the 
form 



fit) « c(r, - ty , 



(8) 



and we study the 'instantaneous' or running estimates 
for 7 and as a function of time. 

The local fits are done as follows: we first produce the 
new function 



din fit) 
dt 



f{t)/f'it). 



(9) 



If fit) is of the form ([s]) then our new function satisfies 
git) K, (T* — t)/"f. Therefore, a linear fit of g(i) will give 
and 7. More explicitly, we have the local expressions 



.f{t)nt) 
fit)' . 



and 



T,(i) = t + 



fit) fit) 
fit) fit) -fit)' 



(10) 



(11) 



The latter local expressions can be used with any suitable 
fit method of the data, not necessarily linear fits. 

In practice, as our time series are given on an 
equally spaced temporal grid, we proceed in the following 
straightforward manner. First we compute ln(/(t)), then 
we use centered finite differences to estimate its deriva- 
tive. Inverting this data furnishes estimates of git) at 
the midpoints. Using again centered finite differences 
produces estimates of I/7 on the original grid, thus al- 
lowing the determination of local estimates for both 
and 7. Note that this algorithm basically amounts to a 
local 3-point nonlinear fit. 

The values of git), T»(i) and 7(t) obtained in this way 
from the TG data for the vorticity supremum ||w||oo are 
displayed in Fig. [sj It is apparent that git) presents 
an inflection point around t = 3.3 corresponding to a 
maximum value of 7 that is above — 1 . Thus local in time 
power-law extrapolations around t — 3.3 are inconsistent 
with the BKM theorem that requires 7 < — 1. However, 
when t is larger than 3.6, the value of 7 goes below — 1 
and thus becomes compatible with BKM. 

On the other hand, there is no sign that the data values 
of 7 and are settling down into constants, correspond- 
ing to a simple power-law behavior. 

Recall (see Section III B ) that the last reliable value of 
||w||oo at resolution 4096"^ is at t = 3.85. Thus, due to 
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FIG. 5: (Color online) Time evolution of a) inverse logarith- 
mic derivative Q at all resolutions (see legend); b) extrapo- 
lated Tlj (111 (solid black line: T, — t) and c) running value 
of 7 (10 1, both only at resolution 4096^ (red circles). 



our 3-point extrapolation method, the last reliable data 
point is at i = 3.825 in Fig. [Sj^a) and at t = 3.8 in Figs. 
[Sj^b) and [sj^c) . The data corresponding to 7 and T* are 
also displayed in Table 



Time 


7 


T, 


3.7 
3.75 
3.8 


-1.42 
-2.06 
-1.04 


4.09 
4.26 
4.02 



TABLE III: Power-law fit parameters 7 and T* (see Eq. ([sjl) 
for the vorticitysupremum ||tt)||cx) determined at resolution 
4096=* (see Fig.jUb) and Fig. ^c)). 

Thus our conclusion for this section is that although 
our late-time reliable data for ||ti;||oo(i) shows ^{t) < — 1 
and is therefore not inconsistent with BKM. clear power- 
law behavior of ||a;||oo(i) is not achieved. 



IV. BRIDGING ANALYTICITY-STRIP 
METHOD AND BKM THEOREM 



< C,^2a+5/2W, Vie [0,T). (12) 
This bound is valid for any e > 0, where 



kezLdUZ?^<,„\{o} 



-3-2 e 



and Qp is defined by 
1 



E |kP^|S(k,OI' 



(13) 



(14) 



Notice that 2r2p is the square of the Sobolev seminorm 

|u(-:*)li/P- 

Motivated by the numerical results of Section [III B| let 
us assume, at a given time t, a behavior of the energy 
spectrum ^ of the type 



E{k) 



(15) 



Notice that n and S are functions of time. When n < 6 
and S tends to zero, this gives a UV-divergence: 



e+5/2 



Plugging this into the bound (12 1, and using the BKM 
theorem, we get J ' 6{t)~^'^T~^dt = 00, where is the 
hypothetical singularity time. 

At this point, again motivated by our numerical re- 
sults, we assume n = const. < 6 and assume a power-law 
behavior for the analyticity-strip width of the form 

S{t) oc (T, - tf , 

where F > is a constant. Replacing this into the above 
integral we conclude that 



(T* -t)(^3+5-^)^di = 



i.e., a finite-time singularity can be attained only if the 
exponents satisfy (—3 + f — e)r < —1 for any e > 0. 
Taking the limit e — > we deduce finally 



A. Motivation and simple estimates 

The vorticity maximum Ht^^l', i)||oo "^^^ found to de- 
crease when the resolution is reduced at any given 
time t > 3.7 (see the above discussion following Fig. 
[ijb)). This strongly suggests that, in this late-time 
regime, ||<^(', t)||oo ^^^^ significant contributions coming 
from high-wavenumbers modes. In this context, the fol- 
lowing short heuristic argument is provided as a motiva- 
tion for the more rigorous mathematical results to follow. 

Consider the well-known Sobolev inequality, which can 
be derived using the same hypotheses as in Lemma [7] 



r > 



6 — rt 

In words: "if the analyticity-strip width S{t) goes 
to zero as a power law, then the exponent must 
be greater than or equal to ." 

The main difficulty to overcome in order to material- 
ize the above heuristic arguments into a firm basis, is 
that the common Sobolev bound ( 12 1 has a problem at 



e = : the constant is equal to infinity there, so tak- 
ing the limit as we did is not fully justified. We provide 
the solution to this problem by finding a new rigorous 
bound, sharper than the common Sobolev bound, and 
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which gives the same optimal exponents without a diver- 
gent constant. 

The second difficulty is that the assumed behavior 
for the energy spectrum (f5l, commonly used in the 



analyticity-strip method, is a very strong condition and 
does not hold uniformly for /c G N. In fact, the evidence 
in analytically solvable models such as the ID Burgers 
equation, is that the behavior ( |15[ ) holds with some expo- 
nents n and S in the region k ^ S~^, (large-fc asymptotic 
limit), and the behavior E{k, t) ^ k~" holds in the region 
1 < k <^ 5^^ , with h < n. We provide the solution to this 
lack of uniformity by introducing a "working hypothesis" 
which is a uniform-in-fc inequality for the energy spec- 
trum, that still retains the spirit of the analyticity-strip 
method. The working hypothesis is verified for the case 
of ID Burgers equation (see below the discussion at the 
end of Sec. 



VII 



B. Mathematical preliminaries 

1. BKM theorem 

We assume the usual hypotheses of the Beale-Kato- 
Majda (BKM) theorem. Let T denote, from here on, a 
generic time so that the velocity field u g C([0, T); H^) n 
C^{[0,T); HP"^) , p > 3, so in particular the quantities 



defined in (14) are bounded for p > 3 : 



np{t) <cp, vt e [o,r) . 

The BKM theorem [Ij states that the assumed regularity 
of the velocity field can be extended up to and including 
the time T if and only if t(T) = ||a;(-, t)||oo dt < oo. 
By 'regular up to and including the time T' we mean 
u e C([0,r];i7P) nCi([0,T];iJP-i), p>3. 

Definition 1. We define the maximal time of regularity 
G (0, oo] as the earliest time for which u ceases to be 
in C{[0,T];HP) n C^{[0,T]; RP-^) , p>3. 
If < oo we speak of a finite-time singularity. 

With this definition, we conclude that the time integral 
appearing in the BKM theorem converges for all T < 
and diverges at T = T*: /J"* ||a;(-,i)||oo dt = oo . 



Definition 2. The kinetic energy spectrum E{k, t) is de- 
fined as the sum of squares of modulus of Fourier coeffi- 
cients over spherical shells 



E{k,t) 



E 

fc-l/2<|k|<fc + l/2 



u(k,i)|^ 



(16) 



It is easy to check that the TG symmetries imply 
that E{Q,t) = E{l,t) = Vi G [0,T,). Numerical 
observations (see reference [TB] and Section III above) 
lead us to formulate the following working hypothesis 
that will be used to bound the energy spectra: 

Hypothesis 3 (Working hypothesis). From here on, we 
will assume that there exist a constant M > and posi- 
tive functions no(i), (5o(i), continuous on [0,r^), such that 
for all times t £ [0, T^) and all k € Tj, k > 2 we have 



E{k,t) < Mfc-"«(*)e-2*=^o(t) 



(17) 



Remarks, (i) The working hypothesis is consistent with 
the hypotheses of the BKM theorem, (ii) The work- 
ing hypothesis is an inequality defined globally in k and 
is not a large-k asymptotic expansion. Furthermore, a 
large- fc asymptotic expansion is typically of the form 
E{k,t) = Ci(i)fc-"i(*) e-2'='^i(*) and has, in contrast to 
( [l7| , a time-dependent constant Ci{t). Nevertheless, 
asymptotic results can be used to establish the working 
hypothesis in special cases such as the ID-inviscid Burg- 
ers equation (see the discussion below, at the end of Sec. 
VI). (iii) The numerically-obtained fits of the analyticity- 
strip method E{k,t) « C(<)fc""(*) e"2fe5(t) ^^^.g similarly 
related to the working hypothesis. Notice that these fits 
are obtained over a finite range of values of wavenum- 
ber fc, so they give only partial information. At early 
times, when the analyticity-strip width S is big so that 
6k ^ 1, one is in the "large-fc asymptotic limit". At late 
times, when S becomes of the order of the highest resolved 
wavenumber fcmax, we have 6k ^ 1 and thus the fits repre- 
sent the "small-fc range". The relations n{t) > no{t) and 
6{t) > 60(1) are required for consistency with the work- 
ing hypothesis. In practice, we will use the numerically 
obtained n{t) and 6{t) to estimate no{t) and 6o{t). 



3. Classification of solutions in terms of regularity 



2. Working hypothesis for energy spectrum 

An implicit assumption of the analyticity-strip method 
is the existence of the Fourier components of the solu- 
tion of the 3D Euler equations. Taylor-Green (TG) sym- 
metries imply that only modes with even-even-even and 
odd-odd-odd wavenumber components are present (see 
Section II B I . The appropriate definition of the energy 



spectrum is thus 



We see from Definition [T] that a finite-time singularity 
is defined by the condition < 00. Combining this with 
the working hypothesis, a finite-time singularity can oc- 
cur only if lim(_>.T^ <5o(0 = 0. Amongst all possible contin- 
uous positive functions 5o(t) that tend to zero as t — > 
we will consider, to simplify the analysis, only power-law 
type of functions. 

Definition 4. A solution of the 3D Euler equations satis- 
fying the working hypothesis (17), is said to have a finite- 



time singularity of power-law type, with power T > 0, iff 
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the working hypothesis admits a function So{t) that be- 
haves, near t = T^., as 

So{t) oc in~tf. 



We saw in the heuristics Section [lV A| that if the energy 
spectrum is of the form E{k,t) « C(t)A:-"'(*) e'^'^''^*), 
then the exponent n{t) must be less than 6 in order for a 
finite-time singularity to occur. This result will be fully 
formalized in Section |IV C[ but first we need to define 
two types of solutions in terms of the behavior of the 
exponent noit) appearing in the working hypothesis. 

Definition 5. A solution of the 3D Euler equations sat- 



isfying the working hypothesis (17 1, is said to be of strong 



regularity iff the working hypothesis admits an exponent 
rio(i) such that liminf 7io(i) > 6. Otherwise, i.e. if all 

the exponents admitted by the working hypothesis satisfy 
liminfrio(i) < 6, the solution is said to be of mild regu- 
larity. 

The reason for the name "strong" is due to the following 



lemma (to be proved in Section IV C I 



Lemma 6. Let a solution of the 3D Euler equations sat- 



isfying the working hypothesis (17), be of strong regular- 
ity. Then the solution has no finite-time singularity. 

This lemma's assertion is basically the same as the 
well-known fact that there cannot be a finite-time loss of 
analytic regularity without loss of regularity . 

This result can be used as a validation test for numeri- 
cal simulations of 3D Euler fluids. If the supremum norm 
of the vorticity is to grow in time without bound, then 
the exponent noit) must be well below the critical value 
6. Fortunately, all reliable numerical simulations that we 
know of pass this elementary test. 



C. Main results linking Beale-Kato-Majda theorem 
and analyticity-strip method 

1. Sharp bound for vorticity 

Lemma 7. Let u(x, t) be a velocity field satisfying the 
Taylor-Green symmetries and with energy spectrum de- 
fined by equation (16). Let u: — \7 x u be its vorticity, 

[0,2 7r]3. Then the 



defined on the periodicity domain D 
following inequality is verified for all times t G [0,r); 



||a;(.,t)||^ < E V^Hk + l)Eik,t)Sk, (18) 



k=2 



where Sk = #{k G Z^^^UZj?^,,, : fc-1/2 < |k| < fc + 1/2} 
is the combined number of lattice points ( of the form odd- 
odd-odd or even- even- even) in a spherical shell of width 
1 and radius k G 



Proof. The vorticity field is defined in terms of its Fourier 
components by w(x, t) = e*'^'^a)(k, t). There- 

fore 



(19) 



for all x G D. The LHS of this equation can be re- 
placed by the supremum norm. Also, we use the identity 
|tD(k,i)| = |k| |u(k,i)| on the RHS and obtain 



|k||S(k,i)|. 



Assuming that u is regular so the above sum over the 
lattice converges, we can rewrite the sum over spherical 
shells of width 1 and radius k G Z+. We get 



oo 
k=2 



I 



\ 



E 



|k||u(k,t)| 



\ **^'=^odd^ even 
\ fc-l/2<|k|<fc-|-l/2 



We proceed to bound the terms in brackets, for a given 
k G First, notice that the highest possible value of |k| 
is equal to \/k{k -I- 1). We obtain the preliminary result 

/ \ 

OO 



\\^{-M, 



< 



Ev/fc(fc+^ 



fc=2 



E 



|u(k,i)| 



\ fc-l/2<|k|<fc + l/2 

Second, the remaining sum in brackets is related to the 
energy spectrum E{k,t), equation (|4|, by virtue of the 
Cauchy-Schwartz inequality. We have 



E |S(k,t)| < V2E{k,t) 



fc-l/2< |k 



<k + l/2 



\ 



fc-l/2<|k|<fe + l/2 



which establishes the Lemma. 



(20) 

□ 



Remarks. The proof is independent of any evolution 
equation that u might satisfy. Only two inequalities have 



been used to get the bound (18 1, and these inequalities 
are quite sharp: 



First, the bound (19) is saturated when all phases are 



equal in the Fourier expansion for the vorticity field at the 
position of vorticity maximum. This saturation indeed 
takes place in one-dimensional systems that blow up in a 
finite time, such as the inviscid Burgers equation (work 
in progress). 



Second, the bound (20 1 is saturated when all the terms 



are equal in the sum over the spherical shell of fixed ra- 
dius k. Physically, such saturation should be observed in 
a fully isotropic scenario, i.e., when the terms |u(k, i)|^ 
depend more on the wavevector's modulus |k| than on its 
direction k/|k|. 

In contrast, the Sobolev bound (12) would be satu- 



rated only for unphysical scenarios where the energy 



10 



spectrum E{k,t) has a compact support in fc-space and 
is independent of the wavenumber k on that support. 



Thus the Sobolev bound ( 12 1 will be less sharp than the 



new bound (18 I for any realistic energy spectrum that 
decays as fc — > oo. 

Practical form. We provide a more practical form of 
the sharp bound (18 1, by noticing that Sk ~ Trfc^ as A: — >■ 
oo. Under the hypotheses of Lemma [7j we readily obtain 
the estimate 



(21) 



fc=2 



where c = 2-^/11/3. This constant was computed by 
direct inspection of the maximum deviation from the 
asymptotic formula Sk ~ tt/c^. Although this estimate 
seems not as sharp as the original one, it will be enough 
for the practical situation where the analyticity-strip 
width d{t) tends to zero and the main contribution comes 
from the 'ultraviolet region' k ^ 1. 



2. Implications of BKM Theorem: General result 

Let us replace the working hypothesis for the energy 
spectrum (17) into the bound (21 1. The sum over k>2 



can be written in terms of the so-called polylogarithm 
function. We obtain the bound 



|la;(,t)IL<cVMLi{'^ 



where Li (s, z) is defined by 



(22) 



(i) Let < z < 1 and let p,q be two non-negative 
numbers. Then we have Li (p, z) < Li (q, z) <;=^ P ^ Q ■ 

(ii) Let < 27r and /ef r G M \ Z_|_. Then 

Li (r, e'') « r(l - r) (-/.)-! + f^fc^/ , 



fe=0 



k\ 



where Q is the Riemann zeta function. 
(Hi) Let \fi\ < 2tt and let s G Z_|_. Then 



,s-l 



Li(s,e^)«^^[i7,_i-ln(-^)]+ , 

k^s-1 



fc! 



where Hp = p-th harmonic number, with 

Ho = 0. 



h=l 



We are now ready to prove 

Lemma [6[ Let a solution of the 3D Euler equations 
satisfying the working hypothesis (17 1, be of strong regu- 



larity. Then the solution has no finite-time singularity. 

Proof. By definition, solutions of strong regularity sat- 
isfy the working hypothesis with liminft_>T, nQ{t) > 6. 
Therefore, using Lemma |9](i) on equation ( p3| , we ob- 
tain /'^* Li (l -I- e,e"''o(*)) = oo , for some e G (0,1). 
Now, using Lemma |9](ii) with ?■ > 1, we obtain that the 
integrand is continuous in time. Therefore = oo. □ 



Li(s,z) = ^fc "z'' = U(s,z) 



k=2 



and Li (s, z) is the Jonquiere's function (or polyloga- 
rithm): Li (s, z) = X]fc=i ^ 



Combining the bound ( 22 ) with the BKM theorem we 
obtain the following 

Theorem 8. Let a solution of the 3D Euler equations 
satisfy the Taylor- Green symmetries and the working hy- 
pothesis |j7p . Then its maximal regularity time T* must 
satisfy 



Li 



(^'^-2,e-^o(*))di = oo. (23) 



Proof. The proof is a direct application of the BKM the- 
orem to inequality ( 22 1 . □ 



At this point it is necessary to state without proof some 
properties of the polylogarithm: 

Lemma 9. The polylogarithm function Li (p, z) satisfies 
the following properties: 



3. Implications of BKM Theorem: Singularity scenarios 

Theorem [8] represents our 'bridge' from analyticity- 
strip method to BKM Theorem: a singularity of the so- 
lution at time can be attained only if the parameters 
no{t) and 6Q{t) satisfy equation (23). 



Recall that for a singularity to occur, the function 
6a{t) must tend to zero as i T*. The polylogarithm 

Li (^^^4^ - 2,e-''o(*)) has a branch point at no(t) = 

6, So{t) = (see Lemma[9](iii)) so the asymptotic behav- 



ior of the integrand ( p3| ) as Sq {t) depends sensitively 
on the behavior of the function ngit) near the 'critical' 
value 6. To avoid this branch point, we introduced solu- 
tions with strong and mild regularity in Definition [5j 

The two following main results exploit the conse- 
quences of Theorem [8] in singularity scenarios. They pro- 
vide us with a criterion on how fast must 6o{t) decay to 
zero in order to sustain a singularity. 

Corollary 10. Let a solution of the 3D Euler equations 
satisfy the Taylor- Green symmetries and the working hy- 
pothesis (17). Let the solution be of mild regularity, i.e.. 
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limmf t^Tt no{t) < 6, where is the maximal regularity 
time. Let limt-j.j'^ 5{){t) = 0. Then, satisfies 



50 (t) 



dt = oo, 



for all in (— oo, liminf no(t)] H (— cxi,6). 



Proof. Let n_ be in (— oo,liminf no(i)] H (—00,6). From 

n_ < liminf(_s.T^ no(i); using Lemma |9] (i) on equation 
( 23 1 we obtain 



Li 



Now, since n_ < 6 and the function 5a{t) tends to 
zero as i — )■ T*, we can use Lemma |9] (ii) to bound 
the integrand Li — 2, e^''"'*)) by a constant times 



, which completes the proof. 



□ 



Finally we consider the hypothetical situation of a 
finite-time singularity of power-law type, as described in 
Definition [i] (5o(t) cx (T, - t)^ , with < cx). 

Corollary 11. Under the hypotheses of Corollary^T^ the 

solution of the 3D Euler equations has a finite-time sin- 
gularity at time < 00, of power-law type with exponent 
F, only if 



F > 



6 — n_ 

for all 7i_ in (—00, liminf no(t)] H {—00, 6). 

Proof. The proof follows directly from Corollary ITO) □ 



V. ANALYSIS OF ANALYTICITY-STRIP 
WIDTH IN TERMS OF BKM THEOREM 



a) 
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FIG. 6: (Color online) Comparison of the bounds for the 
Taylor-Green flow at resolution 4096'^. a) Lin-Log plot, "old 
bound" is the RHS of the inequality ( 12 1, taking e = 0.1 and 
text) 



C, = 3.9 
inequality ( 18 1 



New bound" is the RHS of the sharp 
b) Interpolated time derivative of the loga- 
for a time range localized near the change of 



rithms of a 
trend. Same parameters as in a) 



old and new bounds are not too bad at the beginning of 
the computation (t = 0) , with an initial ratio of 5 : 2 be- 
tween the new bound and the vorticity supremum norm. 
Subsequently, the bounds become increasingly less sharp, 
and the new bound attains a ratio 165 : 1 with the vor- 
ticity supremum norm at t = 4. However the slope of 
the new bound's curve is comparable to the slope of the 
vorticity-supremum-norm curve. 

In order to make a more quantitative comparison of the 
slopes. Fig. [6]jb) shows the logarithmic rates of growth 
for old bound, new bound and vorticity supremum norm. 
In that order, these rates satisfy the ratios 7 : 5 : 4 at the 
resolved time t « 3.85. 



A. Quality of bounds 



B. Analysis of S in terms of BKM 



Several bounds were used in Section |IV| We now pro- 
ceed to test their sharpness, when they are applied to 
the numerical data of Sectio n |III| Fig. |6] shows a com- 
parison of the new inequality ( 18 1, and the old inequality 
(fT2l taking e = 0.1 with = 3.9. Note that he value 



of Ce (see Eq. (13)) can be estimated by the integral 
^7r3-V2e yielding 3.75 at 



V3 



3-2£ 



dk 



.1. A more careful computation of the discrete sum 
gives Ce > 3.9, the value used to generate Fig. |6] 

The data in Fig. [6ja) displays two important facts: (i) 
The new bound is sharper than the old bound through- 
out the computation, particularly at the reliable end of 
the simulation, t ^ 3.7, when the three curves show a 
change of trend and the old bound diverges at a faster 
rate than the new bound (see also Fig. [6|^b)). (ii) Both 



We now proceed to see if the accelerated decay ob- 
served in the decrement S{t) and quantified in Fig. |4|d) 
can correspond to a power-law. To wit, we use the same 
local 3-point method than that described in Section [ill C| 
(see Eqsj9 10 and 11 1. The behavior of g{t) is presented 
in Fig. l7[ and the corresponding T*(t) and T(t) are pre- 
sented in Table |lVl 

The results for exponent and predicted singular time of 
Table [Tv] have to be read carefully. Because of the local 3- 
point method used to derive them from the data in Table 
ITT) they use the values of 5 at t = 3.65, 3.7, 3.75, 3.8, 3.85, 



the last one being marginally reliable (see Sec III B I . In 
fact, they amount to linear 2-point extrapolation of the 
data in Fig. |7](see the inset): is the intersection of the 
straight line extrapolation with the time axis and F is the 
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FIG. 7: (Color online) Temporal evolution of the inverse log- 
arithmic derivative (|9| computed from the same values of S 
as m Fig. ^d); 512^Ybrown triangles), 1024^ (blue squares), 
2048^ ( green diamonds) and 4096'^ (red circles). 
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r. 


T. 




3 ^max 


103 ^max 


3 ^max 


103 ^max 


3.7 


0.283 


0.383 


3.81 


3.83 


3.75 


0.354 


0.393 


3.83 


3.83 


3.8 


1.41 


1.36 


4.00 


3.97 



TABLE IV: Power-law fit parameters F and T, (see Eq. ([8 
for S{t) determined at resolution 4096'^ on full interval 3 



k < fcmax (same as in Fig. J4 
103 < fc < fcmax (see Table In 



and Fig. [7| and on subinterval 



inverse of the slope. One can guess that there is room 
for a power-law type of behavior, with exponent F w 0.4 
if we consider the data a,t t — 3.7, 3.75 and F « 1.4 if we 
include the data a,t t — 3.8. 

We now use Corollary 11 (see Sec. |IV[ ) to test if these 
estimates of power-law are consistent with the hypothesis 
of finite-time singularity. There, the product r(6 — n_)/2 
must be greater than or equal to one if finite-time singu- 
larity is to be expected. With the conservative estimate 
n_ = 3.9 obtained by inspection of Fig. |4|^b) (or equiva- 
lently using the values of n in Table |ll| , we obtain that 
r(6-n_)/2 < 1 for the data at t = 3.7 and t = 3.75, but 
r(6 — n_)/2 > 1 for the data &t t — 3.8. These results are 
insensitive to the fit interval, see Table |TV] Therefore, if 
the latest data is considered. Corollary 11 cannot be used 
to negate the validity of the hypothesis of finite-time sin- 
gularity. However, there is no sign that the data values 
of F and in Table [TV| are settling down into constants, 
corresponding to a simple power-law behavior. 

Another piece of analysis consists of comparing the sin- 
gular time predicted from the data for the decrement d{t) 
with the singular time predicted from the direct data for 
the vorticity supremum norm. They seem both to be 
close to w 4 (compare Table [lV| to Table [mi. 



In this context, we should perhaps mention Feynman's 
rule: "Never trust the data point furthest to the right", 
a comment attributed to Richard Feynman, saying ba- 
sically that he would never trust the last points on an 
experimental graph, because if the people taking data 
could have gone beyond that, they would have. Higher- 
resolution simulations are clearly needed to investigate 
whether the new regime is genuinely a power law and 
not simply a crossover to a faster exponential decay. 

Our conclusion for this section is thus similar to that 



of Sec. HI C although our late-time reliable data for S{t) 
shows F(6 — n_ ) /2 > 1 and is therefore not inconsistent 
with our Corollary 11, clear power-law behavior of d{t) 
is not achieved. 



VI. CONCLUSIONS 

In summary, we presented simulations of the Taylor- 
Green vortex with resolutions up to 4096'^. We used the 
analyticity strip method to analyze the energy spectrum. 
We found that, around t ~ 3.7, a (well-resolved up to 
t ~ 3.85) change of regime is taking place, leading to 
a faster decay of the width of the analyticity strip S{t). 
In the same time-interval, preliminary 3D visualizations 
displayed a collision of vortex sheets. Applying the BKM 
criterium to the growth of the maximum of the vorticity 
on the time-interval 3.7 < t < 3.85 we found that the 
occurrence of a singularity around t ~ 4 was not ruled 
out but that higher-resolution simulations were needed 
to confirm a clear power-law behavior for ||w||oo(t). 

We introduced a new sharp bound for the supremum 
norm of the vorticity in terms of the energy spectrum. 
This bound allowed us to combine the BKM theorem 
with the analyticity-strip method and to show that a 
finite-time blowup can exist only if S(t) vanishes suf- 
ficiently fast. Applying this new test to our highest- 
resolution numerical simulation we found that the be- 
havior of S (t) is not inconsistent with a singularity. How- 
ever, due to the rather short time interval on which S{t) is 
both well-resolved and behaving as a power-law, higher- 
resolution studies are needed to investigate whether the 
new regime is genuinely a power law and not simply a 
crossover to a faster exponential decay. 

Let us finally remark that our formal assumptions of 
Section |IV C| are motivated and to some extent justified 
by the fact that, in systems that are known to lead to 
finite-time singularity, the analogous of the working hy- 
pothesis (17) is verified. For the analogy to apply, a 

This is 



version of the BKM theorem must be available, 
the case of the 1-D inviscid Burgers equation for a real 
scalar field u{x,t) defined on the torus: 



du 



du 
dx 



Vx e [0,27r], e [0,T, 



which admits a BKM-type of theorem [5S], with singu- 
larity time defined by J"^* ||ua;(-, t)||oo dt = oo. 
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In the 1-D case, the analogous of our bound (21) is 



oo 

\\u,i-,t)\\^<dJ2k./E{kJ). 

k=l 

Using the simple trigonometric initial data u{x, 0) — 
sin(a:), the energy spectrum can be expressed in terms 
of Bessel functions that admit simple asymptotic expan- 
sions. It is straightforward to show (see [lOj for details) 
that, for t < = 1, one has the large-fc asymptotic 
expansion 



E{k,t) 



1 



^-3g-25s(t)fe^ 




with 



while, at i = = 1, 

9 fi2/3 

In fact, the k~^/'^ power law appears already before 
(see the remark following Eq. (3-10) of reference |10]). 

It is easy to check that the analytical solution admits, 
for all k and for all t sufficiently close to T*, a working 
hypothesis (IT7|) of the form 



E{k,t) < Mk~'"> exp(-2(5o(i)fc), 

with analytically-obtainable functions no{t) = 8/3 and 
So{t) cx (T, - with r = 3/2. The analogous of Corol- 
lary [TT] gives the inequality 



r > 



4-no ' 



which is saturated by the analytically-obtained expo- 
nents no = 8/3, r = 3/2. 



VII. ACKNOWLEDGEMENTS 



tion IV C Definition |2] and the working hypothesis (Hy- 
pothesis [3j) are modified slightly in the general case. Ac- 
cordingly, the new bounds leading to Lemma [7] and The- 
orem [8] need to be modified slightly to accommodate the 
general case. The crucial derived relations between do 
and riQ in Lemma |6] and Corollaries [lO] and [TT] will ap- 
ply directly to the general periodic case and will not be 
discussed. 

The main technical difference is that the new bounds 
presented in Section |IV C| apply for a flow with TG sym- 
metries (see Section II B I which imply that only modes 
with even-even-even and odd-odd-odd wavenumber com- 
ponents are populated. The general periodic case does 
not follow this restriction, which slightly modifies the 
bounds. We will assume, to simplify matters, that the so- 
called zero-mode of the velocity field is identically zero: 

u(0,i) = 0, ViG[0,r), 



Notice that all remaining wave numbers are populated. 
This means that all sums involving the scalar k in equa- 
tions (18) and (21 1 will start effectively from k = 1. 



Also, because modes with mixed even-odd wavenum- 
ber components are allowed, the definitions of 5*^ in 
Lemma 2 and constant c in equation (21) must be 
replaced by more appropriate quantities. Therefore, 
the corresponding general periodic versions of Lemma [7] 
(equation (18l) and practical bound (equation ( pT] )) are: 



Lemma [7| (general periodic version of Lemma [t]). 

Let u(x, t)be a velocity field with energy spectrum defined 
by equation ^ and let LiJ — V x u be its vorticity, defined 
on the periodicity domain D = [0, 2ttY . Then the follow- 
ing inequality is verified for all times t G [0, T) when 
the sum in the RHS is defined, and independently of any 
evolution equation that u might satisfy: 



M-Moo < E ^j2Hk + l)E{k,t)Sl, (Al) 

k = l 



where 5^ = #{k G : k - 1/2 < \k\ < k + 1/2} is the 
number of lattice points in a spherical shell of width 1 
and radius k G Z_|_. 
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"Singularities in three-dimensional Euler equations: sim- 
ulations and geometry". 



Appendix A: Extension to general periodic flows 

Here we provide the generalization to non TG- 
symmetric periodic flows of the results presented in Sec- 



Practical bound, general case. 

\\u,{;t)\\oo<^'flk^VE{kJ). 
fc=l 



(A2) 



where c' — 6\/2- 

We can easily check that the bounds for Taylor-Green, 



equations (18) and (21), are sharper (by a factor close 



to 2) to their respective general bounds, equations (All 
and ( |A2| . 

Finally, Theorem [8] is replaced by 

Theorem |8j. Let a solution of the 3D Euler equations 
satisfy the working hypothesis |T?|) with k = 1 included. 
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Then the maximal regularity time of the solution must 
satisfy 



10 



[1] J. T. Beale, T. Kato, and A. Majda, Communications in 

Mathematical Physics 94, 61 (1984). 
[2] C. Bardos and E. Titi, Russian Mathematical Surveys 

62, 409 (2007). 
[3] R. Kerr, Physics of Fluids 17, 075103 (2005). 
[4] T. Y. Hon and R. Li, Journal of Nonlinear Science 16, 

639 (2006). 

[5] M. D. Bustamante and R. M. Kerr, Physica D: Nonlinear 

Phenomena 237, 1912 (2008). 
[6] J. Gibbon, Physica D: Nonlinear Phenomena 237, 1894 

(2008). 

[7] P. Constantin, C. Fefferman, and A. J. Majda, Communi- 
cations in Partial Differential Equations 21, 559 (1996). 
[8] J. Deng, T. Y. Hou, and X. Yu, Communications in Par- 
tial Differential Equations 30, 225 (2005). 
[9] J. Deng, T. Y. Hou, and X. Yu, Communications in Par- 
tial Differential Equations 31, 293 (2006). 

[10] C. Sulem, P.-L. Sulem, and H. Frisch, Journal of Com- 
putational Physics 50, 138 (1983). 

[11] C. Bardos and S. Benachour, Annali della Scuola Nor- 
male Superiore di Pisa - Classe di Scienze 4, 647 (1977). 

]12] D. Gottlieb and S. A. Orszag, Numerical Analysis of 
Spectral Methods (SIAM, Philadelphia, 1977). 

]13] U. Frisch, Turbulence: The Legacy of A. N. Kolmogorov 
(Cambridge University Press, 1995). 

]14] U. Frisch, T. Matsumoto, and J. Bee, Journal of Statis- 
tical Physics 113, 761 (2003). 

]15] G. I. Taylor and A. E. Green, Proc. Roy. Soc. A 158, 499 



(1937). 

[16] M. E. Brachet, D. I. Meiron, S. A. Orszag, B. G. Nickel, 
R. H. Morf, and U. Frisch, Journal of Fluid Mechanics 
130 (1983). 

[17] M. E. Brachet, M. Meneguzzi, A. Vincent, H. Politano, 
and P. L. Sulem, Phys. of Fluids A 4, 2845 (1992). 

[18] C. Cichowlas and M.-E. Brachet, Fluid Dynamics Re- 
search 36, 239 (2005). 

[19] E. Lee, M. E. Brachet, A. Bouquet, P. D. Mininni, and 
D. Rosenberg, Phys. Rev. E 78, 066401 (2008). 

[20] A. Pouquet, E. Lee, M. E. Brachet, P. D. Mininni, and 
D. Rosenberg, Geophysical & Astrophysical Fluid Dy- 
namics 104, 115 (2010). 

[21] P. D. Mininni, D. Rosenberg, R. Reddy, and A. Pouquet, 
Parallel Computing 37, 316 (2011), ISSN 0167-8191. 

[22] S. S. Ray, private communication (2012). 

[23] S. S. Ray, U. Frisch, S. Nazarenko, and T. Matsumoto, 
Phys. Rev. E 84, 016301 (2011). 

[24] I. Kukavica and V. C. Vicol, Discrete and Continuous 
Dynamical Systems - Series A 29, 285 (2011). 

[25] M. D. Bustamante, Physica D: Nonlinear Phenomena 
240, 1092 (2011), ISSN 0167-2789. 

[26] This new behavior of the Euler TG vortex is somewhat 
similar to the acceleration in the decrease of 5 that was 
reported in MHD for the so-called IMTG initial data at 
resolution 2048'^ in reference [19] . 

[27] http://www.vapor.ucar.edu 



